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INVESTIGATING  OF  LONGITUDINAL  DEVELOPMENT  PARAMETERS  THROUGH  AIR 


In  this  work  the  simulation  of  the  Extensive  Air  Showers  was  performed  by  investigating  the  longitudinal 
development  parameters  (N  andXmax)  by  using  a  system  for  air  shower  simulation  which  is  called  AIRES  version  2.6.0  at 
the  energy  range  (1014-1019eV)  for  different  primary  particles  like  (gamma,  electron,  positron,  proton  and  iron  nuclei)  and 
different  zenith  angles.  The  comparison  of  simulated  longitudinal  profile  was  fulfilled  for  different  hadronic  models 
(SIBYLL,  QGSJET99  and  SIBYLL  S16). 
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INTRODUCTION 

When  high-energy  cosmic  ray  particles  penetrate  the  Earth  atmosphere  they  interact  and  generate  a  cascade  of 
secondary  particles,  that  is  called  Extensive  Air  Showers  (EAS),  but  the  only  way  to  get  information  on  these  particles  is  to 
study  EAS  cascades  which  they  initiate  [1].  When  a  cosmic  ray  enters  the  atmosphere  it  collides  with  the  nucleus  of  an  air 
atom,  producing  a  number  of  secondary  particles.  These  go  on  to  make  further  collisions,  and  the  number  of  particle 
grows.  In  the  end,  the  energy  of  the  shower  particles  is  cut  down  to  the  point  where  ionization  losses  to  overcome  and  their 
number  start  to  descend  [2]. 

The  longitudinal  development  of  EAS  is  depended  on  the  energy  and  type  of  the  incident  primary  particle  [3]. 
One  of  its  characteristics,  the  atmospheric  depth  of  shower  maximum  (Xmax),  which  is  often  used  to  reconstruct  the 
elemental  composition  of  primary  cosmic  rays.  The  number  of  charged  particles  (N)  in  EAS  as  a  function  of  atmospheric 
depth  is  intimately  related  to  the  type  and  energy  of  the  primary  particle  which  can  be  simulated  by  AIRES  code  [4]  .A 
realistic  air  shower  simulation  system  which  includes  electromagnetic  interaction  algorithms  and  links  to  the  hadronic 
models  SIBYLL  [5]  and  QGSJET  [6]. 

In  this  work  the  investigation  of  the  longitudinal  development  profiles  in  EAS  was  performed  by  using  AIRES 
code  for  different  hadronic  models  like  (SIBYLL  [5],  QGSJET99  [6]and  SIBYLL  S1.6[7])  packages  for  the  simulation  of 
hadronic  processes.  The  simulation  was  performed  for  different  primary  particles  (gamma,  electron,  positron,  proton  and 
iron  nuclei)  at  the  high  energy  range  (1014-1019eV)  for  different  zenith  angles  (0°,  10°and  30°). 

Longitudinal  Profile 

The  Gaisser-Hillas  formula  gives  the  approximate  number  of  charged  particles  as  a  function  of  atmospheric  depth 
along  the  shower  axis  [8] : 


SHOWER  SIMULATION  BY  DIFFERENT  HADRONIC  MODELS 


AL-RUBAIEE  A.  AHMED  &  AHMED  JUMAAH 


Al-Mustansiriyah  University,  College  of  Science,  Department  of  Physics,  Baghdad,  Iraq 


ABSTRACT 


(1) 


136 


Al-Rubaiee  A.  Ahmed  &  Ahmed  Jumaah 


where  Xis  the  atmospheric  depth  (in  g  cm"  );  Nmax  is  the  number  of  charged  particles;  Xmax  is  the  depth  of  shower 
maximum  and  X  is  a  characteristic  length  parameter  (a  scale  constant  with  a  value  70  g  cm"2)  [9].  This  is  a  Gamma  function 
form  which  naturally  arises  in  cascade  theory,  and  assumes  that  the  first  interactionis  at  X=0.The  first  interaction  point  is 
the  location  of  initial  collision  of  the  cosmic  ray  particle  with  atmosphere.  A  forth  parameter  is  often  introduced  into 
Eq.(l),  ostensibly  to  allow  for  a  variable  first  interaction  point  [10]: 


(x  —  Xr,  v  iXwm-xky* 
 -M  exp  (X™*  -  X)/X. 


(2) 


Where  X0is  the  depth  of  the  point  of  first  interaction.  The  value  ofX0  depends  on  the  collision  cross  section  and 
hence  the  energy  and  mass  composition  of  the  particle.  The  depth  of  shower  maximum  Xmax  depends  on  the  position  of  X0, 
the  shower  energy  and  composition.  But  one  can  express  the  shower  longitudinal  development  as  function  of  the  shower 
agesthat  defined  as  s=3X/X+2Xmax  instead  of  depth  X,  by  using  s.  Translating  the  depth  X„mx  into  age  s  and  using  the 
normalized  shower  size  n=N/Nmax,  theEq.  (2)  becomes[ll]: 


n(s)  =  (l 
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whereTm=Xmax  IX  and  T0=  X0  /X. 
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In  this  work  we  apply  the  AIRES  simulation  code  to  study  generation  of  primary  particles  like  (gamma,  electron, 
positron,  proton  and  iron  nuclei)  induced  air  showers  at  primary  energy  (1014  to  10I9eV)  and  explores  the  longitudinal 
development.  Although  the  simulation  of  air  showers  at  the  lowest  energy  is  not  a  very  realistic  application  we  intend  in 
the  present  work  to  investigate  the  longitudinal  development  by  the  different  hadronic  models(SIBYLL,  QGSJET99 
and  SIBYLL  S16)  at  the  highest  energy  ranges  (1014-1019eV)  for  different  zenith  angles.  All  simulations  performed  by  the 
AIRES  code  were  obtained  by  using  the  thinning  factor  10"4relative  thinning  case. 

Figure  (1)  shows  the  longitudinal  development  for  the  primary  particles  (gamma,  electron  and  proton),  when  air 
shower  curves  represent  number  of  particles  with  fixed  atmospheric  depth  (in  g  cm"2).  The  results  of  simulation  are 
compared  between  different  hadronic  models  at  the  energy  range  (1014-1019eV)  for  vertical  showers. 
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Figure  1:  The  Longitudinal  Development  Simulated  by  AIRES  Code  for  Vertical  Shower  and  Different  Hadronic 
Model  (SIBYLL,  QGSJET99  and  SIBYLL  S16),  for  Different  Primary  Particles 

The  longitudinal  development  for  the  primary  particles  (electron,  positron  and  iron  nuclei)  which  displayed  in 
figure  (2)  reflects  another  importance  or  characteristic  of  the  EAS  development  that  simulated  by  AIRES  code  for  #=10°, 
and  the  fluctuations  are  very  large  for  10"4  relative  thinning  level,  that  they  reduce  immediately  when  the  thinning  is 
lowered. 
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Figure  2:  The  Longitudinal  Development  Simulated  by  AIRES  Code  for  0=10°  and  Different  Hadronic  Models 
(SIBYLL,  QGSJET99  and  SIBYLL  S16),  for  Different  Primary  Particles 


138 


Al-Rubaiee  A.  Ahmed  &  Ahmed  Jumaah 


The  figure  (3)  Shows  the  results  of  simulated  longitudinal  development  in  EAS  for  different  primary  particles 
(proton,  gamma,  iron  nuclei),  at  the  energy  range  (1014-1019eV),  for  (9=30°.  Through  this  fi  gure  one  can  see  the  effect  of 
the  hadronic  models  on  the  fluctuations  of  the  longitudinal  development  of  primary  particles,  in  the  same  conditions  as  in 
figure  (2)  for  different  hadronic  models. 
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Figure  3:  The  Longitudinal  Development  Simulated  by  AIRES  Code  for  0=30°  and  Different  Hadronic  Models 
(SIBYLL,  QGSJET99  and  SIBYLL  S16),  for  Different  Primary  Particles 

CONCLUSIONS 

In  the  present  work  the  code  AIRES  was  used  for  simulating  the  development  of  air  showers  in  the  atmosphere 
initiated  by  different  primary  particles  (gamma,  electron,  positron,  proton  and  iron  nuclei)  at  the  energy  range 
(1014-1019eV)  with  different  zenith  angles.  The  simulation  of  the  longitudinal  development  has  shown  an  opportunity  of 
primary  particle  identification  and  knowledge  of  the  atmospheric  depth  for  the  cosmic  ray  spectrum. 

The  obtained  results  led  to  the  conclusions  that:  (a)  Found  difference  between  calculations  in  different  hadronic 
models;  (b)  Showers  induced  by  heavy  particles  are  more  sensitive  to  atmospheric  depth  maximum.  The  main  advantage  of 
the  given  approach  consists  possibility  of  a  longitudinal  development  analysis  of  real  events  which  detected  with  extensive 
air  shower  array  and  reconstruction  of  the  primary  cosmic  ray  energy  spectrum  and  mass  composition. 
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